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^ . Abstract 
•*-> . 

Let f2 be an operator semigroup with generator A in a sequentially complete locally convex 
topological vector space E. For a semigroup with generator A + D, where D is a bounded linear 
operator on E, two integral equations are derived. A theorem on continuous dependence of a 
semigroup on its generator is proved. An application to random walk on Z is given. 
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00 : 

£NJ ■ Let ft, ft' be operator semigroups [H [2j [3] (in this article we consider only one-parameter semi- 

\q . groups) with generators A, A' respectively and let D = A' — A. It is known [3, p. 77], that under 

O ' appropriate assumptions 

fl' = ft + flD * ft 1 , (1) 
ft 1 = ft + ft' * Dfl, (2) 

where * signifies convolution, in this case - of operator-valued functions. We call these equalities the 
perturbation formulae. 

In the present paper we prove the perturbation formulae under more general than in [3] assump- 
tions, then deduce from them the strongly continuous dependence of a semigroup on its generator 
and, finally, apply these results to the approximate calculation of operator semigroups in coordinate 
spaces. The idea of the use of the perturbation formulae is to regard the equalities ([I]) and ([2]) as 
equations for the unknown function ft'. 

An important example of a semigroup is the transition probability of Markov process. For this case 
the perturbation formulae were discovered and are systematically exploited by Portenko [H [5] . One 
can find ibidem also a generalization of these formulae for an inhomogeneous process (its transition 
probability creates a two-parameter semigroup). 

Denote the space where the semigroups are defined by E. Below, A\ stands for A' and Aq stands 
for A (likewise ft\ and fto). We impose the usual in the theory of operator semigroups assumptions: 

(i) E is a sequentially complete locally convex topological vector space; 
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(ii) Ai (i G {0, 1}) is densely denned linear operator such that for some Ao G M the resolvent 
i^(A ) = (A 1 - Ai)' 1 is defined on the whole E and the family {(Ai^(A)) w : A > l,m G N}, is 
equicontinuous. 

We remind that the equicontinuity of a family {/q, 6 G 0} of maps from E to E means [2] that for 
any continuous on E seminorm | • | there exists a continuous seminorm || • ||' such that for all G G 
and q G E \\fg(q)\\ < \\q\\' . 

Theorem 1. Let conditions ^ and (jn]) be fulfilled and the operator D = A\ — Aq be defined and 
continuous on the whole E. Then for any t > 

n 1 (t) = fi (i)+ I fi (t-r) J DOi(r)dr, (3) 
J o 

n 1 (t) = n (t)+ [ fti(t-T)L>ft (r)dT. (4) 
Jo 

Proof. Since D is defined everywhere, Aq and A\ have the same domain which we denote S. 

Let, at first, condition (Jn]) hold for some Ao < 0. Then the Hille-Yosida theorem [2j IX. 7] asserts 
that, firstly, the operators Ai generate equicontinuous semigroups 0^ strongly continuously depending 
on t and, secondly, Ri(X) = J °° e~ xt Qi(t) dt. Consequently, equality ([3]) holds for all t > iff the 
equality 

i?i(A) = R (X) + R (X)DR 1 (X) (5) 

holds for all A > Ao- To prove the latter we note that by the definition of resolvent Aj-Rj(A) = 
Ai?, ; (A) - 1 = Ri{X)Ai on 5 so that R DR l = RqA^ - RqAqRx = i? (Ai?i - 1) - (XR - 1)R X . This 
proves for the restriction of its both sides on S. Assumption (Ju|) allows to extend it to E. 

In case Ao > the same argument applied to the operators Ai = Ai — X$l yields (J3J) for f2j thereafter 
it remains to note that, obviously, f2j(i) = e _A °*J7j(t). 

Swapping A\ and Aq, we convert © to dH). □ 

As the first application of the perturbation formulae we will prove that the correspondence between 
a generator and the generated semigroup is continuous w.r.t. the strong topology. The distinction of 
this result from the Trotter-Kato theorem ([21 IX. 12], [61 IX. 2. 16]) is that the condition for convergence 
of a sequence of semigroups is formulated in terms of generators rather than resolvents. 

Theorem 2. Assume that: the space E satisfies condition and the operators A n (n G Z+) satisfy 
condition (jn]); for each n G Z + the operator B n = A n — A is defined and bounded on the whole E; for 
any q G E 

B n q -» (6) 

as n — > oo. Then for arbitrary t > 0, q G E and any seminorm \\ ■ \\ among those generating together 
the topology of the space E 

lim sup ||f2 n (s)g — f2(s)g|| = 0. 

s<t 

Proof. From ([6]) we deduce by the Banach-Steinhaus theorem [7], Theorem 2.6] that 

sup \\B n \\ < oo. (7) 

n 
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Denote g n = £l n q, w n = ||n n ||, W n (t) = sup s < t w n (s) (similarly g, w, W), V n = \\B n \\W n . Substi- 
tuting in ([2|) Q' by f2 n and D by B n , we get by the properties of seminorm 

w n (s) < w(s) + V n (s) / w n (r) dr, (8) 



o 



su V \\g n (s) - g(s)\\ <W n (t) [ \\B n g(r)\\ dr. (9) 

s<t JO 

Having iterated inequality (jHJ) once, we get 

w n (s) < w(s) + V n (s) J (w(r) +V n (r) J w n (p)dp\ dr 

r>S 

\2 



<w(s)+V n ( S ) w(T)dT + V n (s) 2 (S - P )w n (p) dp. 

Jo Jo 

Iterating further, we obtain in the limit 

w n (s) < w(s) + V n (s) f e v ^ s){s - T) w{T) dr, 

whence 



o 



W n (t) < W(t)e 



\B n \\W(t) 



The last equality together with ([7|) shows that sup n W n (t) < oo. Now we deduce the conclusion of the 
theorem from and (|6|), relying on the dominated convergence theorem applicable due to ([7|). □ 

The general scheme of applying Theorems 1 and 2 is as follows. Assume that the semigroup with 
generator Ao can be calculated explicitly (in infinite-dimensional spaces this situation is rare but 
possible — a nontrivial example will be given below) and there exists a strongly converging to the 
generator A sequence (A n ,n G N) of generators such that for each n the operators B n = A n — A and 
D n = A n — Ao (and therefore A — Ao) are defined on the whole space and continuous. If herein the 
operators D n are sufficiently simple (in some sense simpler than A — Ao), then equation ([3]) or (j4j) 
can prove quite manageable, which will be illustrated below. It enables us to explicitly calculate the 
semigroup £l n for any n. Theorem 2 asserts that, for sufficiently large n, we thus obtain a satisfactory 
approximation of the semigroup £1. Unfortunately, this theorem (as well as the above-mentioned result 
of Trotter and Kato) does not give the rate of convergence and, consequently, the exact rule for the 
choice of n. So, it is expedient, while applying the theorem, to calculate some initial segment of the 
sequence rather than a single semigroup. In this case it is simpler to perturb on the nth step 
j4 n _i rather than A, that is to put D n = A n — A n _\. 

Let us refine the obtained results for vector topological spaces of numeral sequences on a finite 
or countable set X (coordinate spaces). Let E be some coordinate space and E' be its dual. We 
denote likewise a linear continuous operator A on E and its adjoint operator on E', discerning them 
by the argument position: to the right from A in the first case and to the left in the second one. If E' 
consists not only of numeral sequences, but the subspace E' = E' n is invariant w.r.t. the action 
of the adjoint operator, then we consider the latter as an operator on this, coordinate restriction of 
the dual space. Then both the operator A and its adjoint are determined by the function on X 2 (their 
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coordinate representation): 

{Aq){x) =J2 A (x,y)q(y), q G E, 

y 

(pA)(y) = ^2p(x)A(x,y), p G E' . 

X 

In the terminology of differential equations the semigroup SI generated by the operator A is the 
fundamental operator, or matrizant, of the system of differential equations for numeral functions 
q(x, •): 

q(x, t) = J2 A (x, y)q(y, t), x £ X, (10) 

y 

which is tantamount to the single equation 

q = Aq (11) 

for an IR^-valued function. 

Let us find the coordinate representation of the fundamental operator in the case when E = l\ (X) 
(the choice of the space will be explained below) and A is bounded, that is, 

ll^ii = su P J2\ A ( x >y)\ < °°- ( 12 ) 

* X 

Obviously, the solution of equation (lllh with initial condition q(0) = qo € h is an Zi-valued function 

and 

< lko||ie l|A|1 *, (13) 

where || • ||i is the norm in the space l\. 

Let G(-, y, ■) denote the solution of system (11Q|) with depending on a parameter y initial condition 

G(x,y,0) = 5(x,y) (14) 

(8 is Kronecker's function). According to (fTT|) and (JT3J) sup^ Yl x \G(x,y,t)\ < \\ A\\e^ A ^ . This together 
with (fl"2j) means that for any qo € l\ the series ^2 y G(x,y,t)qo(y) admits term-by-term differentiation. 
So we obtain the following equation for its sum q(x,t): 

q(x, t) = Y^ G(x, y, t)q (y) =J2Y1 A ( x > V> 
y v i 

Permuting the summations, we become convinced that the functions q(x, •) satisfy the system of 
equations (fTUl) . This jointly with the obvious equality q(x, 0) = qo(x) shows that G(-,-,t) is the 
coordinate representation of the operator £l(t). We denoted it otherwise than the operator itself 
because one can regard the system (|10[) as a single differential-difference equation for q(-, •), in which 
case G is naturally interpreted as Green's function. 

A similar argument shows that G{x, •, •) is the solution of the conjugate to (|10p system 

p(y,t) = ^2p(x,t)A(x,y), y £ X, (15) 

X 

with depending on a parameter x initial condition flHJ). 
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The choice of coordinate space is determined by mathematical or physical reasons which may 
be irrelevant to functional analysis. For example, in probability theory, the differential equations 
for semigroups Q = AQ, £1 = 0,A correspond to backward and forward Kolmogorov's equations for 
the transition probability matrix of a homogeneous Markov chain. In this case it is naturally to set 
E = /^(X), where X is the state space of the chain. Then the dual E' has l\ as its subspace and 
E'q = l\. This choice of the space together with additional assumptions about A allows to interpret 
the solution of (|15p as the one-dimensional distribution of the chain. In statistical physics, the set 
(fT5|) is called (again under some assumptions about A) a master equation and p has the physical 
meaning of concentration (density). The latter must be bounded but need not be summable. In this 
case E should be matched in such a way that E' = loo- This is ensured by the choice E = l\ (so that 
E' Q = E'). In what follows we consider just these spaces. 

Denoting 

F(x,y) = J2 G(x, r))D(r),y), (16) 
v 

H(x,y) = ^2D(x,0G(Z,y), (17) 

e 

we rewrite equalities (pQ) and ([2j) in the coordinate form 

G'(x, y) = G(x, y) + ^2 F & * G"(£, v), (18) 

G'(x, y) = G(x, y) + ^2 G'(x, r,) * H(rj, y). (19) 
v 

These are countable sets of integral equations with respect to the functions G(x, y) of argument t 
suppressed in notation. In the general case they are not simpler than (jlOp and (|15[) . But if the operator 
D or its adjoint is finite-dimensional then in order to solve them it suffices to find the resolvent of a finite 
system of integral equations. Indeed, denote Hi = {77 : 3£ D(£,rj) / 0}, H2 = {£ : 3r/ £)(£, 77) 7^ 0}, 
so that F(x,rf) = as 77 ^ Hi, H(£,y) = as £ ^ H2. Restrict in (PT8j) the range of x by Hi and in 
(fTUj) the range of y by H2. Then for every y (fT8|) is a system of equations w.r.t. the functions G'(x, y), 
x € Hi, of argument t; for every x (|19p is a system of equations w.r.t. G'(x,y), y € H^- The kernels 
of the equations and, consequently, the resolvents do not depend on y in the first case and x in the 
second one. Having found one of the two resolvents, we express via it those functions that enter the 
system, thus converting equality (fT8|) or (fT9f) to the explicit formula for the rest of G'(x,y). 

The idea of the use of the perturbation formulae is that, according to Theorem 2, the solution of 
an infinite system of differential equations (jlOp with bounded A can be approximated with a desired 
accuracy by the solution of its finite subsystem. Even if the operator A is unbounded, but we are able 
to solve the set of equations 

q (x, t) = ^2 Mx, y)qo(y, t), x G X, (20) 

y 

with such an operator Aq that the operator A — Aq is bounded and both A and Aq satisfy condition 
([n]) , then the solution of (|10p can be approximated with a desired accuracy by the solution of the set of 
equation finitely perturbed from (|20p . The accuracy can be evaluated by the perturbation formulae. 
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As an example we consider a spatially inhomogeneous random walk on Z. Forward Kolmogorov's 
equation for the probability p(y, t) for the walking particle to be at state y at time t has the form (|15|) 
with operator A which has the following nonzero entries: 

A(x, x + 1) = A(x), A(x, x — 1) = /i(x), A (a;, x) = — A(x) — yu(x), 

where A(x) and /u(x) are the intensities of transitions from x to x + 1 and x — 1, respectively (we 
assume that the functions A and \i are bounded). For the operator A corresponding to spatially 
homogeneous walk (A(x) = fJ,(x) = 1) the transition probability is given by the well-known formula 
[8, Formula II. 7. 7] 

G {x,y,t) = e-%_ J/ |(2t), 

where I n is the modified Bessel function of order n. Thus, for any random walk whose transition 
intensities are distinct from unity only for a finite set of states, we are able to calculate the transition 
probabilities exactly. To this end we pass to Laplace transforms and solve in the Laplace domain 
the set of equations (|18j) which under the restriction x € H will be a finite set of linear equations. 
The same formula (|18p will give the answer in Laplace domain for any x, y € Z. The inverse Laplace 
transform of the result, that is the function G'(x,y,t), turns out to be a finite linear combination 
of modified Bessel functions of argument 2t whose orders do not exceed \x — y\ (the coefficients are 
algorithmically computable) . 
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